On an exact hydrodynamic solution for the elliptic flow 
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Looking for the underlying hydrodynamic mechanisms determining the elliptic flow we show that 
for an expanding relativistic perfect fluid the transverse flow may derive from a solvable hydrody- 
namic potential, if the entropy is transversally conserved and the corresponding expansion "quasi- 
stationary", that is mainly governed by the temperature cooling. Exact solutions for the velocity 
flow coefBcients V2 and the temperature dependence of the spatial and momentum anisotropy are 
obtained and shown to be in agreement with the elliptic flow features of heavy-ion collisions. 
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CL|! I. INTRODUCTION 
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' The hydrodynamic description of the formation and development of a quark-gluon plasma (QGP) in high-energy 
' heavy-ion collisions has met a considerable success 1]. In particular, the hydrodynamic features seem to be, at least 
partly, required in order to take into account the second Fourier coefficient V2 of the transverse flow of particles, the 
so-called elliptic flow. One writes 0, [3] for the azimuthal multiplicity distribution 

^ = ^{l + 2v2Cosi2^) + ---}, (1) 
1—^ dip Zn 

O ; 

^ ' discarding for simplicity other Fourier coefficients non- relevant here. The experimentally observed values of V2, which 
are due to the anisotropy of the initial state collisions at nonzero impact parameter, are sizeable enough to require 
important collective effects of particle production. These are better reproduced by hydrodynamical properties of the 
^sj I flow in some early stage of the quark-gluon plasma formation. 

J> ' The theoretical estimates of the elliptic flow are obtained from numerical studies based on various versions of the 
, hydrodynamic models. Indeed, a full study requires not only to deal with the solution of the relativistic hydrodynamic 
\[ • equations but also with the definition of appropriate initial conditions and a model for the mutation of the QGP pieces 
of fluid into particles. The numerical studies (c/. reveal that the QGP as a fluid is "almost perfect" since its 
viscosity is remarkably weak, even if the model dependence may account for some variation on the quantitative 
estimates. This observation has a considerable theoretical impact, since it points to a strongly coupled plasma, 
guiding a large part of theoretical interest towards strongly coupled gauge field theories. 

Our goal in the present work is to try and identify by explicit analytic solutions the basic hydrodynamic mechanisms 
at work in the elliptic flow. For this sake, we have to simplify (or even idealize) the description of the QGP formation 
in a heavy-ion collision while keeping the main physical ingredients. Among other simplifications that we will discuss 
now, our study will assume the QGP to be a perfect fluid without viscosity. We will also restrict our analysis to the 
; ^ ■ transverse flow in the central rapidity region where the hydrodynamic description is better suited. 
Ci I The main characteristic feature of the hydrodynamic description of QGP formation in heavy-ion collisions appears to 
be a nontrivial combination of: a) the large longitudinal momentum and energy boosts provided to the created medium 
by the initial state, b) the (presumably fast) equilibration of the energy density and all three pressure components 
due to local thermalization required by hydrodynamics. As a matter of fact, the first stage of the hydrodynamical 
description of particle production in high energy collisions is mainly governed by the longitudinal fiow, that is the 
expansion of the relativistic fluid in (1-1-1) dimensions. Already the pioneering papers of the hydrodynamic approach 
d, 1H| based their analysis on this property. 

In the mean time, the 4-dimensional hydrodynamical feature of the system is kept with the thermodynamic relations 
which, through the local temperature T and the Equation of State (EoS), lead to 
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where r is the proper-time, Cg the speed-of-sound, which we will assumed to be constant in the following, and e,p the 
energy and (isotropic) pressure density, respectively. Our remaining notations are s for the entropy density, a;^={o,- -3} 
for the space-time coordinates {xq = t), and u^-{o,- -3} for the fluid 4- velocity (with lower indices) in the Minkowski 
metric t]^'^ , with signature (1,-1,— 1,-1) satisfying the normalization condition 

w^u^ = Uq - U3 — u^^ = 1 ; = uf + u\ . (3) 

Our idea for studying the transverse motion of the fluid is that it is also driven by the longitudinal evolution, but 
slowly enough to be considered as "quasi-stationary" , that is in such a way that its time evolution is essentially 
related to the temperature cooling. Indeed, the seed of transverse momenta is indirect and there should be, at least 
during some first stage of the hydrodynamic evolution, no strong back-reaction on the longitudinal motion. The 
"quasi-stationarity" hypothesis will allow for an exact solution for the elliptic flow and in general for the hydrodynamic 
regime in the transverse plane. To be concrete we state the following conjectured properties of the hydrodynamic flow: 

• a) Tranversally isentropic. Since the overall entropy should be conserved, we will conjecture that the transverse 
flow is itself (approximately) isentropic, i.e. we write the following equation 

[dxAsUl) + da=AsU2)]transverse = « ' (4) 

• b) Quasi — stationary : Since the time dependence of the transverse entropy distribution is absent from ([4]), we 
will close the equations for the transverse flow by assuming that its hydrodynamic evolution is smooth enough 
to be driven only by the temperature change. We thus consider the equation and solutions of a temperature- 
dependent stationary flow, with a source emitting a fluid at a given temperature (and thus transverse speed, 
see further). 

Hence, in the regime when the longitudinal expansion is dominant, the transverse motion is conjectured to be 
smoothly driven by the overall local temperature of the fluid, which provides the 4-dimensional feature^ of the system 
through the thermodynamic relations 

The plan of the paper is the following: in the next section II, using the hydrodynamic potential for a stationary flow 
, we derive the analytic equation obeyed by the azimuthal distribution of entropy and thus the elliptic flow. Then, 
in section III, we find the exact solutions for the transverse flow. In section IV we apply our solution showing that 
the obtained elliptic flow retain good qualitative features observed in reality or in realistic numerical hydrodynamic 
model studies. A discussion of the hypotheses and our conclusions and outlook form the final section V. 

II. AZIMUTHAL ENTROPY DISTRIBUTION 

As we shall see now, the conditions a) and b) lead to nontrivial properties of the fluid and to analytic solutions for 
the elliptic flow. It is obtained from an hydrodynamic (KK) potential [7] derived by Khalatnikov and Kamenshchik 
for a stationary transverse isentropic flow. The "quasi-stationary" hypothesis allows us to extend its applicability to 
a slow transverse motion of the fluid and to find the general exact solution of the elliptic flow. In fact, the existence 
of a hydrodynamical potential obeying a linear equation is known since long Q for the longitudinal evolution. 
Recently Iff], it was possible to express interesting analytic solutions for the entropy distribution dS/dy where the 
(hydrodynamic) rapidity is deflned by y = ^ log(uo-l-U3)/ log(uo — U3). We shall follow the same method as in for 
the transverse flow case, and find the general solution of the KK potential in order to obtain the azimuthal distribution 
of the entropy dS/d(p giving access to the elliptic flow. 

However, one crucial difference of the transverse w.r.t. the longitudinal case is the velocity-temperature relation 
between uq, the time component of the velocity (and thus also u± — ^Mq — 1, the modulus of the transverse one) to 
the local temperature. This comes from the relativistic Bernoulli relation, verified by a stationary fluid [ll|, namely 

Tuo = T^l + ul = To . (5) 



In that respect, our picture is different from a purely transverse hydrodynamic flow Q. 
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In the same conditions, the whole evolution from some initial temperature to the final freeze-out one is constrained 
to be in the supersonic regime as we will verify later through our equations. The condition writes 



v± 



1/2 



> c.. 



Hence, the isentropic transverse evolution starts at a given temperature Tj such that 

Ti<Ts = To v/l-c2, 



(6) 



(7) 



and the velocity increases when the temperature decreases from T/, reaching eventually ultra-relativistic values before 
hadronization. For convenience, we will from now on introduce the variable 



/ = i log 



ilog 



log v±. 



(8) 



The derivation of the KK potential comes briefly as follows [7] : Together with the transverse entropy conservation 
the equations for the transverse flow close with the projection to the transverse plane of the energy- momentum 
conservation relation 9^T„„n = 0, again by neglecting the time derivatives w.r.t. the transverse gradients. After 
nontrivial transforms, presented in appendix \A\ one obtains the system of equations 

dxiisui) + dxAsu2) = 

d^,{Tu2) ~ dx,{Tu{) - 0. (9) 

Then, using the "hodograph" 0, H, H, E^j [HI inversion of variables {xi,X2) {lif) and combining Eqs.Q and ([9]), 
one arrives at the formulae expressing the kinematic (now dynamical) variables {xi, X2) in terms of the hydrodynamic 
variables through a suitably defined KK potential function xif, I), namely 




a{ip, I) = If + arctan 





m 






dx 




dip 





(10) 



where we have parameterized 



Ui 
Xi 



U± COS ip 

x\ cos a 



U2 = u±_ sm Lp 
X2 = x±_ sin a 



(11) 



The KK potential function x{Vj is solution of a linear equation obtained by closing the hydrodynamic equations 
system using the EoS 



1 



1 



dl 



(12) 



Note the zero coefficient at e' = Cs which signals the supersonic bound ([61 [7|) at T = Ts . In fact the system expands 
in the vacuum for T < Tg, while it is compressed when T > Ts, cf. Q . Hence the physical solutions are restricted to 
the supersonic range T < Tj <Ts. 

The KK potential and its equation have been reproduced from Q. The calculation of the entropy distribution is 
now parallel to the one [lO] (see _12,] for an early version) used in the (1-1-1) dimensional case. Considering an entropy 
flux normal to the tangential line element {dxi, dx2), one has to compute 



dS — su2dxi — suidx2 ■ 

Using the formulae pH)) for the expression of the line element in terms of the potential, one has 



(13) 



6 d\ 

- — = xi COS + X2 sin (f = x± cos(a — ip) 



To dl 
e~' dx 



X2 COS <p — x\ sin <~p — x^ sin(a — ip) , 



(14) 
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which, by differentiation with respect to I and 0, gives 



dS ■ 



T 



dldip d(p 



dl 



d'^X ^ dx 
dip^ dl 



dip 



(15) 



where we used the relation u±e~'' = uq = Tq/T. 

At fixed temperature (and thus fixed I), which is the case considered further on, one gets the azimuthal entropy 
distribution 



dS 

dip 



T 



dl 



sT 



To(l-e27c2) 



dx{v,l) d^xifJ) 



dl 



dP 



(16) 



where the second expression comes from the KK potential equation (H^)- Note again the singular coefRcient at Tg, 
corresponding to the lower bound of temperature. 



III. EXACT SOLUTION OF THE TRANSVERSE FLOW 



In our idealized hydrodynamic framework, without hadronization, one relates the entropy distribution to multiplic- 
ity, dS/S ^ dN/N. Hence, the elliptic flow is defined by the azimuthal entropy distribution (fTB]) through a Fourier 
expansion similar to ([1]), namely 



dS_ 

dip 



and thus 



1. 

2^ 



V2 



{1 + 2i;2 cos{2ip) 



/dy>cos(2y;)f (^) 



} 



(17) 



(18) 



The eccentricity can be obtained as a function of temperature (or I) through Eqs. (fT4|) in terms of the KK potential 
xi'^J) as: 



(a^i - xl) ^ J dip{xl{ip>,l)-xl{ip,l)) 
~ {xl + xD ~ J dip{xl{ip,l)+xl{ip,l)) 



Jdip< 


cos2(^ 


'(Ox' 




+ 2sin2^ lltj 


J dip 







(19) 



Note the characteristic feature of the hodograph method: a geometrical parameter, here e, is expressed in terms of 
dynamical ones, here the temperature. Once the solution is found, one has to invert these relations in order to restore 
the hierarchy between "cause" and "effect" . 

We can now proceed by looking for the general solution resulting from the KK potential solution of (fT2|) . For this 
sake, it is convenient to expand the potential 



X(¥', = /3o(0 + E /^pW cos(2p^) 



in Fourier coefficients Pp{l) which verify the equations 



(e^' - 1) /3;'(0 + e'\c;' - l)/3;(0 - (c;^ e^' - l) (3,(1) = 



(20) 



(21) 



where primes denote derivatives with respect to I. 

As it is well known, the solution is in general a suitable combination, with constant coefhcients, of two independent 
solutions of the second-order equations (|^. Using standard textbooks, one finds for p ^ 



.(1) 



^ ^p ^2. 2 



,21 



5-c: 



-2 2 



5-c: 



Cs — 1 



-2 2 



(22) 
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while we single out the first component cq Po{l), whose derivative simply writes 



In (|22p . 2-^1(0, b; c, z) is the usual hypergeometric function while G2 2 ^ 



,2i 



(23) 

denotes a Meijer function 13]. 



The function (S'^^ {resp. P^^) is the regular {resp. irregular) solution^ at e' = of (|2ip . Hence, any arbitrary 
combination of these two independent functions is a solution of (j2ip . The boundary conditions will define the specific 
linear combinations to be chosen for the general solution which may be obtained from the Green function of the 
problem. For illustration, in the special case p = 0, one writes 



+ CXD 



e(-z) l-e 



Foil) di, 



(24) 



where -fo(0 describes a distribution of sources in temperature convoluted with the Green function, which in this case 
is just 6(— /3o(0 from ((23| . A straightforward but more tedious expression can be written for all values of p but is 
skipped here for brevity. A specific realization for the sake of our physical problem will be discussed in the application 
section. 

Inserting the general solution (I20l22p for the KK potential in Eq. for the azimuthal entropy distribution, one 
finds 



where we denote 



p, and 



mi) 

= A, and 



(25) 



■(T) = 



2/3^(0 [2m+m)]+j:7=i [Mp+mmiP+i)ii)+'^pmPp+iii) - 2(p+ i)/3p+i(0/3;(o -/3;(0/3;+i(0] 



(26) 

for the eccentricity. It is useful for further use to note that, in the "elliptic approximation" i.e. when one stops the 
Fourier expansion (|20p at p = 1, the eccentricity can be expressed using the same functions and parameters with V2 
as in Eq. (P5)l . namely 



siT) 



2/3^,(0 [2/3i(0 + /3U0] _ 


Wo ' 


2(M^^+A/3f)) 


+ /3«'+A/3f)'' 




2(3',\l) + [4(3fil)+f3?il)]~^ ' 






)V + aM^)') 


2" 



(27) 



Thanks to the analytic solutions we obtained, all the expressions contain an explicit dependence in temperature. One 
should only specify which temperature is physically relevant, e.g., Tj for the initial spatial eccentricity and some 
freeze-out temperature Tf for the observed V2. Using our formulae, one may discuss the dynamical hydrodynamical 
process through the temperature dependence of both the spatial and momentum average anisotropy of the lump 
of quark-gluon plasma. For this sake, we note an interesting parameter-independent relation between the spatial 
eccentricity at any temperature T/ and V2 at any temperature Tf, namely 



2/3i(T/) 



e{Ti) Pi{TiW{Tf) 1 + 



2/3i(T,) 



1-2£2(T,) 



1 j_ 

\^ ^ 2l3i(Ti) J 



(28) 



For p = 0, Eq. II21I I is only first-order for /3Q(i) and thus introduces only one arbitrary coefficient cq. 
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A general physical comment is in order about the parameters A and p defining the relevant solutions in the "elliptic 
approximation" . Using a source of given temperature, the parameter A, which corresponds to the relative strength 
of the two independent solutions of the second-order differential equation ((2T|) . will specify the initial condition of 
the elliptic flow. The parameter p, which is geometrical in nature since it gives the relative strength of the elliptic 
harmonic in (|20p . will be related to the initial centrality of the reaction. For general initial conditions, the more 
general Green function formalism, c/. (|24|) . has to be used. 



IV. EXACT ELLIPTIC FLOW: APPLICATIONS 



Taking into account the linear equations for the potential and entropy distributions, c/. (|12|15p . the determination 
of the elliptic flow boils down to defining properly the boundary conditions, i.e. the sources of the hydrodynamic 
expansion, which are given functions of temperature and azimuth. In the following we will assume that the source is 
simply given by a delta-function at the initial temperature T/ of the process and a given initial eccentricity profile. 
We fix it by the condition that V2{Tj) = while s{Tj) is maximal. Note that the solution satisfies the constraint 
T/ < Ts, i.e. the fluid is always supersonic. 

In the "elliptic approximation" for which the Fourier expansion of the potential (|20p is limited to the two first 

orders, the observables V2 (|25p and e (|27p depend only on two relevant parameters, namely p ~ c^^''/co, obtained from 
the Fourier expansion (j20p and A = cp''/cj"'^\ that is the coefficient ratio between the regular and irregular solutions 
of the potential equations (|12l2ip . 



Determination of X. From the previous discussion, A is chosen in such a way that V2{Tj) = 0, where s{Tj) is maximal. 
As an illustration of the discussion, the temperature dependence of both V2 and s that we obtain with our definition 
of the initial condition, is displayed in Fig[Tl for a given value of the geometrical anisotropy parameter p ~ 0.8. The 
value of Ti is lower but nearby the speed-of-sound lower limit of temperature Tg . 
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FIG. 1: (Color online) Compared temperature dependence for the momentum V2{T) and the spatial e{T) anisotropies. The 
curves correspond to the initial temperature source ai T = Tj (see text). The dashed line is for the supersonic lower bound Ts- 
The geometrical anisotropy parameter (see text) is p = 0.8 and the speed of sound is the reference one Cs = l/VS. 



Determination of p. The determination of the geometrical parameter p, the first anisotropy coefficient of the potential 
X, see is governed by the centrality. In Fig[5]we display e{p) which shows a quasi-linear behavior. This is in good 
agreement with the observed feature of the experimentally reconstructed eccentricity with an observed proportionality 
relation with the centrality c ~ Npart/Nmax, where Npart is the number of participant nucleons. Indeed, one expects 
a simple relation between p and c , up to a rescaling p/ pmax ^ Npart /Nmax- With such a choice and using Ea. (j25p . 



7 



0.8n 




0.0 0.2 



0.4 0.6 0.8 
P 



FIG. 2: e as a function of the geometrical anisotropy parameter p. The observed dependence qualitatively reproduces simulations 
of e as a function of centrality c ~ Npart/Nmax (see text). In this figure we used = 1/3 and pmax ~ 0.8. 



one finds a simple proportionaUty rule of V2 with centrality, namely 



V2 = Pmaxi'^ 



(29) 



which is also expected from hydrodynamical simulations [T4|. Note that, in this framework, pmax is indeed independent 
of the evolving temperature ratio T/Tq, but it may depend on the initial conditions such as the type of heavy- ion 
reaction and the initial c.o.m. energy (or Tq). The T/Tq dependence, given in ([^^ by the function within brackets, 
is uniquely defined from l|25p . In our calculations, the linearity of the formula (|25p for V2 in terms of the normalized 
second Fourier coefficient /? is a direct consequence of the equation for the KK potential which, together with the 
azimuthal entropy distribution, is diagonalized by the Fourier expansion (j20p . It is clear that the formulation of the 
initial eccentricity profile depends on the initial conditions, and we take the curve in Fig[2] as an example. 

Fig[T]is interesting also from the point of view of the dynamics of elliptic flow. Indeed, it is known from hydrodynamic 
models [T^ that the momentum anisotropy, represented in our "quasi-stationary" approximation by the temperature- 
dependent V2{T), rapidly increases as a function of proper time, and thus with decreasing temperature, to reach its 
observed value. It is therefore confirmed to be a good indicator of the early stage of the hydrodynamic expansion. On 
the same footing, the spatial anisotropy, represented by s{T), decreases as the system expands, even reaching negative 
values, i.e. changing the sign of the spatial anisotropy. We observe, in Fig[l] that the transversally isentropic flow 
follows the same qualitative path as a function of temperature cooling. It is also interesting to note that the flnal 
value of V2 (and thus the value of V2/£, where e = e(T/) is the initial eccentricity) is reached rather early and rather 
independently of the choice of the initial temperature for the transverse flow. 

In order to restore the time variable through its dependence on the temperature, we shall make use of a convenient 
rescaling of the temperature equivalent to the expansion time, similar to the one proposed in p^ . where the ratio 
V2(t — to)/£(to) with initial time tq is displayed for different values of impact-parameter and various values of the 
speed-of-sound Cg. One makes the rescaling substitution 

i^-^o) ; ^ = 1/7^ + 7^ ■ (30) 




R 

where R gives an appropriate average estimate of the expanding size of the plasma. In our temperature-dependent 
scheme, we deflne an analogous rescaling using the thermodynamical relation ^ by choosing a "rescaled time" variable 
defined in terms of the temperature as 




(31) 



In Fig[3]we show the theoretical results for V2/S as a function on the rescaled variable 6 for the solution we considered. 
Let us comment both sides of the figure. In the upper graph, the figure displays the dependence on centrality via 
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FIG. 3: (Color online) v^je as a function of the "rescaled time" 9. Up: dependence on the centrality via the p-parameter at 
fixed Cs = l/\/3; Down: dependence on the EoS via cl = (1/3, 1/5, 1/7, 1/10) at fixed eccentricity e — 0.8. For this sake, one 
is led to choose, respectively, p = (0.8, 1, 2, 4), by tuning the values of pmax- The "reduced time" is defined by Eg. pTjl . 



P — Pmaxi^ — c). The value of Pmax has been chosen fixing V2{t — To)/e(To) to match with some realistic value (see, 
e.g. [il, [ij]). One observes the general trend of the 9 evolution as a function of increasing centrality (or decreasing p). 
This trend, which has been empirically observed in hydrodynamic simulations [Tsj is here explained by the nonlinear 
e-dependent correction to V2/e (see (pS)) ). It is easy to realize that the remnant e-dependence in is such that 
it increases for increasing e (the denominator is smaller) and thus it decreases with centrality. Hence our scheme 
reproduces, at least qualitatively, a trend as a function of impact-parameter observed in hydrodynamical simulations. 

Also, in the lower graph of Fig. [3] we display the speed-of-sound dependence of V2{9)/e. Thanks to the "time" 
rescaling (j30p . it is possible to superimpose the different curves, provided an adequate choice of p ensures the constant 
initial s{Ti). As also seen numerically, the analytical dependence over Cg of our resulting formula ([28]) gives a decreasing 
value of the ratio V2{0)/e{Tj) with decreasing speed-of-sound. Here also one finds the observed behavior [15]. However, 
this hierarchy is obtained at rather larger 6 than observed in [T5| . We will comment on that feature in the next section. 

V. CONCLUSION, DISCUSSION AND OUTLOOK 

Let us briefly summarize our results: using the conjecture of a "quasi-stationary" and transversally isentropic 
hydrodynamic regime governing the transverse flow, and for a given EoS, we arrive at a closed system of hydrodynamic 
equations which can be solved by a hodograph transform xi^2 —^T,ip. Thanks to the potential method J| can 
formulate the general solution and give explicit analytic expressions for the hydrodynamic features of the transverse 
flow. In an application to a source with given temperature and constant effective speed-of-sound Cs, we are able to 
give a complete analytic solution. The applications to the determination of the features of the elliptic flow are in good 
qualitative agreement with the observed (or numerical) characteristics: the temperature dependence of the spatial and 
velocity anisotropy (see Fig[T|), the linear behavior of V2 with centrality (c/. Eq. (P5)) ) for a realistic initial eccentricity 
(see Figl21) and the centrahty and speed-of-sound dependence of the ratio V2/£ (see Figl3]). 

Now, possessing an analytic solution, it is possible to go back to the initial assumptions and discuss their range of 
validity. In other terms we may address the question to which approximation can we consider our closed system of 
transverse equations to be a good approximation of the full hydrodynamic equations. To quantify this approximation 
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a meaningful comparison is to give estimates of two quantities which are relevant for the discussion of the two 
hypotheses, a) of a "transversally isentropic" and b) a "quasi-stationary transverse" flow. 

In order to test our conjecture a) and looking to ([4]), we are led to consider the following ratio of entropy flow 
gradients 



dr{suo) drisua) drriT) ~ V 9t(smo) 



(32) 



where V = average, temperature-dependent, expansion rate. Indeed, this ratio governs the effect of the 

time-gradient compared with a typical transverse one. Note, however that the overall transverse entropy gradient is 
zero, by virtue of ([4]). 

In ([5^ . we have replaced the kinematical variables by their temperature-dependent averages defined by our solution. 
The approximation range of a "transversally isentropic" is thus related to the value of the (analytically known) 
expression (j32p to be larger than 1 in a significant range of "reduced time" . In the upper graph of FigU] one sees that 
the transverse over longitudinal entropy gradient becomes indeed significantly larger than 1 for high enough speed of 
sound. For small speed of sound this requires longer "reduced time" . This could explain the features of FiglH low, 
with a "retarded" ordering w.r.t. [l5l |. 

In order to test the consistency of the "quasi-stationary" approach b), in the lower graph of Figdiwe present the 

expansion rate itself V = ^§^^^fy: where the functions t{T), x±{T) are analytically obtained from their definition 
within our temperature-dependent scheme, namely from ([2]) and (fTO|) respectively. Note that this rate is also appearing 
in the denominator of , which shows that the two hypotheses of a "transversally isentropic" and "quasi-stationary 
transverse" fiow are indeed connected, since a slow motion gives rise to a high transverse over time typical entropy 
gradient. From Figl4] we see that both the transversally isentropic and quasi-stationary hypotheses are consistent 



c'=1/3 




FIG. 4: (Color online) Comparison of entropy and kinematic gradients. The analytically known quantities dx^{su±) /drisuo) 
and V = dTix±)/dTT) (see text) are plotted as a function of the reduced time 0. 



at not too short "reduced times" and not too small speed-of-sound. Thus, these hypotheses give a qualitative 
analytic understanding of the transverse fiow. Our qualitative picture seems consistent. However, the time gradient 
is not negligible w.r.t. the transverse derivative, indicating, at least within the initial conditions we choose, that a 
quantitative agreement could be more difficult to be obtained. 
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Another topic is the range of vahdity of our approximation in transverse space. Indeed, due to ([6]), the quasi- 
stationary approach is only vahd in the supersonic dilatation regime, which requires a large enough transverse velocity 
v± > Cs- This could limit the range of validity of the hydrodynamical flow which has been observed only at small 
transverse momentum. It could also compromise the dominance of the longitudinal Bjorken flow determining the 
thermodynamical relations ([2]). We think that this limitation, which should be taken into account for a quantitative 
study difficult to perform analytically, will not endanger the qualitative but explicit solution we found. The study 
of the implications for the transverse momentum dependence of the elliptic flow deserves per se a study which goes 
beyond the scope of the present work, where no mass relation between fluid velocity and transverse momentum has 
been introduced. 

As an outlook, it will be interesting to develop the study of hydrodynamical mechanisms generating the elliptic flow 
by the investigation of other phenomenological aspects, such as the abovementioncd p± dependence, the effect of a 
weak viscosity, etc.... In order to reach more quantitative features, it will be useful to refine the definition of the initial 
conditions. In fact, it could be worthwhile to typically defining a priori the dependence of e as a function of p or 
centrality, and finding the corresponding initial conditions'^ by inverting e.g. ([24|) . In particular, to examine whether 
they could identify more definitely the hydrodynamical mechanisms. On a more theoretical ground, the existence 
of rather simple mechanisms may facilitate the search for a relation to the fundamental gauge field theory, and in 
particular the gauge/gravity dual approach of the elliptic flow. 

All in all, our results suggest an analytic approach to the transverse motion of the fluid, which can clarify the 
behavior of the elliptic flow obtained from the data or numerical simulations. This is related to an approximate 
"quasi-stationary" and "transversally isentropic" property, of the transverse flow for which the time dependence of 
the system comes mainly through the temperature. 
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APPENDIX A: QUASI-STATIONARY TRANSVERSE FLOW OF A PERFECT FLUID 



In this section we derive the basic equations determining the quasi-stationary transverse flow of a perfect fluid. The 
energy-momentum tensor of a perfect fluid is 

T^'" ^ {e + p)u''u'' ^ pr]"" , (Al) 

where e is the energy density, p is the pressure and (/i = {0, 1, 2, 3}) is the 4- velocity in the Minkowski metric r/^", 
with signature (1, —1, —1, —1). It obeys the equations 

5^TJ' = ^ u,d^[ie+p)u''] + ie+p)u^d^u,-d,p = 0, (A2) 

with u^u'^ = 1, and thus 

u^df.u" = 0. (A3) 

From now on and for simplicity, denotes d^. Multiplying the equations of motion (jA2p by Uj^, i.e projecting them 
on the direction of the 4- velocity, and using (jXsp we acquire: 

a^[(e+pK]-M''a^p = 0. (A4) 

Finally, re-inserting (|A4|) into (|A2|) we obtain: 

(e -I- p)u^d^Ui, = -UyU^df_,p + d^p. (A5) 



^ We thank C.Gombeaud for this suggestion. 



11 



Relation (|A5|) holds for a general perfect fluid. For a stationary flow it gives rise to Bernoulli equation HDl, namely 

TuQ = To. (A6) 

Let us now focus on the stationary transverse flow 7], namely setting U3 = 0, which is the regime of interest of 
the present work. In this case the equation of motion (|A5[) for gives d^p — and thus the various quantities do 
not depend on the longitudinal coordinate 0:3. Therefore, e, p and the velocities are functions of xi, X2 only. The 
equations of motion (|A5P boil down to: 

di[{e+p)ul+p\+d2[{e+p)uiU2] = 
di[{e+p)uiU2]+ d2[{e+p)ul+p\ = 

(9i [(e +p)uiMo] + ^2 [(e +p)-U2Uo] = 0, (A7) 



with Uq = 1+1*^ + 1*2. Equations (|A7p can be expressed in terms of the temperature and the entropy density. 
Considering vanishing chemical potential we have: 

p + e = Ts ; de = Tds ; dp = sdT. (A8) 

Using relations (jA8[) . equations (|A7[) become: 

di [{Ts)ul] + sdiT + 82 [{Ts)uiU2] = (A9) 
5i [{Ts)uxU2] + 82 [(Ts)u2] + sd2T = (AlO) 
di[{Ts)u^uo\+d2[{Ts)u2Uo\ = 0. (All) 

Now, the Bernoulli equation (|A6p allows to transform (jAlip to the entropy conservation in the transverse plane: 

5i(sui) + 92(5^2) -0. (A12) 

Finally, equations (|A9p and (jAlOp . with the use of (IA3p . (|A6p . (jAlip . give rise to the same equation: 

di{Tu2)^d2{Tui). (A13) 
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